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Abstract. Recently, the first example of two-neutron decay from the ground state of 
an unbound nucleus, *®Be, was seen (A. Spyrou, et. ah, Phys. Rev. Lett. 108 102501 
(2012)). Three-body methods are ideal for exactly treating the degrees of freedom im¬ 
portant for these decays. Using a basis expansion over hyperspherical harmonics and 
the hyperspherical R-matrix method, we construct a realistic model of '®Be in order to 
investigate its decay mode and the role of the two-neutron interaction. The neutron-'^'Be 
interaction is constrained using shell model predictions. We obtain a ground state for 
*®Be that is under-bound by approximately 0.7 MeV with a width of approximately 0.17 
MeV. Eor such a system, an attractive three-body force must be included to reproduce the 
experimental ground state energy. 


1 Introduction 

Exotic nuclei can be found across the nuclear chart, from light two-proton and two-neutron halo 
nuclei, such as '^Ne and ®He, to heavy two-proton emitters, such as "^^Fe, '^^Ni, and ^'*Zn. These nuclei 
can be modeled as three-body systems, a core plus two neutrons or protons, and three-body wave 
functions can be used to investigate their properties as well as to give insights into their decay modes. 
(See for example |[Tl.) The difference between dinucleon decay (where two correlated neutrons or 
protons decay from a nucleus) and three-body decay (where two uncorrelated nucleons decay from a 
nucleus) is of particular interest. 

Two-proton decay was first theorized in 1960 and was then first experimentally verified over 
40 years later in experiments such as a, a, and 0. These systems have been analyzed in terms of 
three-body models, where knowledge of their decay mechanism has led to a better understanding of 
their structure |[T|. Although there has been much progress in the field of two-proton decay, much less 
is known about two-neutron radioactivity la. Recently, two-neutron decay from the continuum has 
been investigated theoretically I?!. 

An experiment performed in 2012 at the National Superconducting Cyclotron Laboratory first 
measured two-neutron radioactivity from the ground state of the unbound nucleus '®Be 0. Although 
'®Be is an ideal candidate for two-neutron radioactivity (the lowest experimentally measured state in 
'^Be is energetically inaccessible to one-neutron decay), there was controversy over this result when 
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it was first published Q. In the original work |[ 8 |, the authors used two extreme models to explain the 
'®Be decay - one without a neutron-neutron interaction to model three-body decay ['®Be ^ *^Be + n 
-I- n] and one with two strongly correlated neutrons to model dineutron decay ['®Be ^ '"^Be + (2n)]. 
However, in using a three-body formulation, we can exactly treat the degrees of freedom relevant 
for this decay - better understanding two-neutron decay and the neutron-neutron correlation in the 
nucleus. 


2 Solutions to the three-body scattering problem 

To solve the Schrodinger Equation in a three-body model, interactions between the three bodies must 
be defined. Typically, experimentally measured levels in the core -H n system (here, '^Be) would be 
used to constrain the interaction between these two bodies. However, only one level in *^Be has been 
measured, at 1.8 MeV fTOl . and therefore, shell model calculations are used to constrain this inter¬ 
action. Shell model calculations for '^Be im give a d 5/2 state at 2.8 MeV, an ii /2 state at 4.0 MeV, 
and a 6 / 3/2 state at 7.0 MeV. In order to match the experimentally measured state in '^Be, these shell 
model levels were lowered by 1 MeV in the current work. In the three-body model, the depth of three 
Woods-Saxon potentials (one for each of the s,p, and d partial waves) is adjusted to reproduce the 
shell model levels. 

The neutron-neutron interaction is taken to be the GPT interaction M- Although this interaction 
is simpler than others such as AVI 8 d and Reid soft-core d, it reproduces NN scattering ob¬ 
servables up to 300 MeV, which is more than sufficient for the energy region considered in this work. 
The GPT interaction has often been used in the study of two-neutron halo nuclei (see for example, 
d, d, ini, and d). A three-body interaction can also be included in order to adjust the final 
resonance energy to match the experimental value. 

For the three-body system, hyperspherical coordinates, p and 9, are used instead of the standard 
Jacobi coordinates HD- This coordinate system allows the three-body wave function to be separated 
into hyperradial and hyperangular parts (where the hyperangular function can be constructed through 
a linear combination of known functions): 


x^JM 


i J,ri^(p)J/i"(G5,o-i,o-2,^). 

P Ky 


( 1 ) 


This expansion introduces the hypermomentum, K, related to the relative angular momenta, and ly 
Hoi , //max determines the size of the model space. 

Then, all that must be solved is the system of coupled hyperradial equations 
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with scattering boundary conditions 


AxKyUP) + Yj ^KK'rY'(p)xi'y'ip) = 0 , 

' jfK' 


( 2 ) 


x'^y'^^‘{Kp) ^ Kp) - KP)] , (3) 

where k - . The index y contains the quantum numbers / (total orbital angular momentum of the 

two neutrons relative to the core), S (total spin of the two neutrons), j (total angular momentum of the 
two neutrons relative to the core), I (the spin of the core), (relative orbital angular momentum of the 
two neutrons), and ly (relative orbital angular momentum of the core -H (2n) system). The subscript i 
labels the incoming channel. 
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Figure 1. Phase shifts as a function of 
three-body energy of the K = Q channel for 
the largest values of used in the 
calculation. 


Equation Q is solved using the hyperspherical R-matrix method lfT9l . Through the R-matrix 
method, phase shifts can be calculated, allowing for the extraction of resonance energies and widths. 

3 Preliminary Phase Shift Analysis and Future Work 

In this work, energies are extracted from the phase shifts in the channel /T = 0 at the points where 
the energy derivative of the phase shift has a maximum. In Figure we show the convergence 
of the these phase shifts. This system required a large model space for convergence, including 80 
hyperradial basis states, hypermomentum greater than 64, and over 80 hyperangular basis states. 
For values greater than /Tmax = 64, numerical inaccuracies were induced due to the strong repulsive 
centrifugal barrier. By hnding the peak of the derivative of the phase shifts in Figure [1] the resonance 
energy was calculated to be 2.07 MeV. 

This ground state energy found above is under-bound by about 0.7 MeV when compared to the 
experimental ground state energy of 1.35 MeV ID. Traditionally, three-body forces in three-body 
models account for the additional binding coming from degrees of freedom that are excluded from the 
model space. Here, an attractive three-body force must be added to the model in order to reproduce 
the experimental ground state of '®Be. 

A rough analysis of the K - 0 phase shift shows that the width is approximately 0.17 MeV. 
Due to the coupled-channel nature of these calculations, the interpretation of the phase shifts is not 
straightforward. Because of effects from the couplings between the channels, phase shifts do not rise 
sharply through 90° at the resonance energy. Therefore, it is not clear whether the ground state of 
the calculated system is given simply by the lowest resonance energy extracted from a single phase 
shift (as computed here) or if these phase shifts have to be weighted and combined, in the case of 
several resonance energies around the same region. One such method for this weighting would be to 
calculate the widths of the state through an integral relationship, such as that used in ll20l . A study 
comparing the various approaches is underway. 
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